Large Scale Structure Constraints for a Class of f{R) Theories of Gravity 
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Over the last few years much attention has been given to the study of modified gravity theories in 
order to find a more natural explanation for the late time acceleration of the Universe. Nevertheless, 
a comparison of the matter power spectrum predictions made by these theories with available data 
has not yet been subjected to a detailed analysis. In the context of f{R) theories of gravity we 
study the predicted power spectra using both a dynamical systems approach for the background 
and solving for the matter perturbations without using the quasi-static approximation, comparing 
the theoretical results with the SDSS data. The importance of studying the first order perturbed 
equations by assuming the correct background evolution and the relevance of the initial conditions 
are also stressed. We determine the statistical significance in relation to the observational data and 
demonstrate their conflict with existing observations. 

PACS numbers: 04.50. Kd, 04.25. Nx, 95.36.+X 



I. INTRODUCTION 

Despite more than a century of careful scrutiny. Gen- 
eral Relativity (GR) remains the best fundamental the- 
ory of physics for describing gravitational interactions. 
When applied to cosmology, assuming that the large- 
scale geometry of the Universe is described by the 
Friedmann-Laimatre-Robertson- Walker (FLRW) metric 
and adopting a fluid description for the matter content 
consisting of baryons, Cold Dark Matter (GDM) and ra- 
diation, GR gives rise to a set of field equations which 
can be solved exactly to give the simplest expanding Uni- 
verse model - the so called Friedmann model, governing 
the dynamics of the cosmological scale- factor a{t). This 
model has been remarkably successful, giving for example 
the correct abundances of the lightest elements and ex- 
plaining the origin of the Cosmic Microwave Background 
Radiation (CMBR). In the last two decades, however, ad- 
vances in observational cosmology appear to suggest that 
if one wishes to retain the FLRW metric, the Universe 
must have undergone two periods of accelerated expan- 
sion. The first period is needed to explain the flatness 
problem and the near-scale invariant spectrum of temper- 
ature fluctuations observed in the CMBR, while the sec- 
ond period explains the dimming of distant type la super- 
novae relative to Einstein-de Sitter Universe model. In 
order to explain these periods of acceleration, the strong 
energy condition (/i + 3p > 0) needs to be violated. In 
the case of inflation, this is achieved by introducing a dy- 
namical scalar field, while the present day acceleration is 
most easily explained with the introduction of a positive 
cosmological constant. The resulting description of the 
Universe has become known as the ACDM or Concor- 
dance Model [1]. Although this beautifully simple phe- 
nomenological model appears to fit all currently available 



observations (Supernovae la [2], CMBR anisotropics [3], 
Large Scale Structure formation [4], baryon acoustic os- 
cillations [5] and weak lensing [6]), it is affected by signif- 
icant fine-tuning problems related to the vacuum energy 
scale and this has led in recent years to efforts to explore 
alternatives to this description of the Universe. 

The leading alternative to the ACDM model is based 
on modifications of the standard Einstein- Hilbert action. 
This is due to the fact that these changes naturally admit 
a phase of late time accelerated expansion (an early uni- 
verse inflationary phase is also possible [7]). In this way 
Dark Energy can be thought of as having a geometrical 
origin, rather than being due to the vacuum energy or 
additional scalar fields which are added by hand to the 
energy momentum tensor (see [8] and references therein 
for an extensive presentation of the state of the art of 
this programme of investigation) . 

One of the simplest extensions of GR is based on grav- 
itational actions which are non-linear in the Ricci cur- 
vature R and/or contain terms involving combinations 
of derivatives of R [8-14]. An important feature of these 
theories is that the field equations can be written in a way 
which makes it easy to compare with GR. This is done by 
moving all the higher-order corrections to the curvature 
onto the RHS of the field equations and defining an "ef- 
fective" source term, known as the curvature fluid. Once 
this has been done the strong energy condition can be 
easily violated and this gives rise to a curvature fluid- 
driven period of late-time acceleration. Unfortunately 
this comes at the cost of having to study a considerably 
more complex set of field equations, making it difficult to 
obtain both exact and numerical solutions which can be 
compared with observations. Many studies of the expan- 
sion history of f{R) gravity have been performed using 
a range of strategies for numerically integrating the cos- 



mological equations and these studies have highhghtcd, 
among other things, how sensitive results are to initial 
conditions, the presence of rip singularities and oscilla- 
tions in the deceleration and snap parameters, the ex- 
istence and stability of the Einstein Static models and 
bounce solutions [15]. 

Theses difficulties can be reduced somewhat by using 
the theory of dynamical systems [16], which, with careful 
choice of dynamical variables, provides a relatively simple 
method for obtaining exact solutions (via the equilibrium 
points of the system) and a description of the global dy- 
namics of these models for a given f{R) theory [17]. 

Another useful (but more limited) approach is to as- 
sume that the expansion history of the Universe is known 
exactly, and to invert the field equations to deduce what 
class of f{R) theories give rise to this particular cos- 
mological evolution [18]. This has been done recently 
for exact power-law solutions for the scale factor, corre- 
sponding to phases of cosmic evolution when the energy 
density is dominated by a perfect fluid. It was found 
that such expansion histories only exist for modifications 
of the type i?" [19]. It was also shown in [20] that the 
only f{R) theory of gravity that admits an exact ACDM 
expansion history is standard General Relativity with a 
positive cosmological constant and if one wants to ob- 
tain this behaviour of the scale factor for more general 
functions of R, additional degrees of freedom need to be 
added to the matter sector. A more extensive analysis 
of reconstruction methods has been carried out in [21] 
to obtain theories which give an approximate description 
of deceleration-acceleration transitions in cosmology and 
also in [22] where the reconstruction method was based 
on standard cosmic parameters instead of specifying the 
time evolution of the scale factor. 

Because it is possible to find background expansion 
histories which are consistent with the standard ACDM 
model, it is necessary to investigate the growth of struc- 
ture in order to break this degeneracy. This requires 
extending the standard theory of cosmological perturba- 
tions for GR to f{R) gravity. This has been done using 
both the metric based approach [23-26] originally devel- 
oped for GR by Bardeen [27] and the 1 + 3 covariant 
approach first introduced by Ellis and Bruni [28]. For 
example, in [29], evolution equations were obtained for 
scalar and tensor perturbations for f{R) gravity and ap- 
plied to the spatially flat, matter dominated solution of 
i?" gravity given by a{t) = aoi^<^+™' (which is a saddle 
point G of the corresponding dynamical system for these 
theories) . 

The results obtained demonstrate that the evolution 
of scalar perturbations is determined by a fourth order 
differential equation rather than a second order one, so 
that the evolution of density fluctuations contains, in 
general, four modes rather than two and can give rise 
to a more complex evolution than what is obtained in 
General Relativity. It was also found that the perturba- 



tions depend on the scale for any value of the equation 
of state parameter of standard matter (while in GR the 
evolution of the dust perturbations is not scale depen- 
dent) and that there is a characteristic scale-dependent 
signature in the matter power spectrum [30] . It was also 
found that the growth of large density fluctuations can 
occur also in backgrounds in which the expansion rate 
is increasing in time. This surprising result is strikingly 
different from what one finds in GR and could provide 
a strong constraint on some f{R) theories using the in- 
tegrated Sachs- Wolfe (ISW) effect and the matter power 
spectrum [31]. 

These features can be interpreted by comparing the 
system of fourth order equations, which produced them, 
with those governing scalar perturbations for two inter- 
acting fluids in GR [32], because they have the same 
structure, i.e., there are friction and source terms due 
to the interaction of the two effective fluids. On very 
large and on very small scales, the system of equations 
become independent of fc, so that the evolution of the 
perturbations does not change as a function of scale and 
the power spectrum is consequently scale invariant. On 
intermediate scales the interaction between the two flu- 
ids is maximised and the curvature fluid acts as a rela- 
tivistic component whose pressure is responsible for the 
oscillations and the dissipation of the small scale pertur- 
bations in the same way in which the photons operate 
in a baryon-photon system. This suggests that the vari- 
ables describing the fluctuations in the curvature fluid 
can be interpreted as representing the modes associated 
with the additional scalar degree of freedom typical of 
/(i?)-gravity. In this sense the spectrum can be explained 
physically as a consequence of the interaction between 
these scalar modes and standard matter and result in a 
considerable loss of power for a relatively small variation 
of the parameter n. 

The results described above were obtained assuming 
the simple background model G and can therefore only 
provide hints about what features to expect in the mat- 
ter power spectrum when realistic expansion histories are 
considered. In order to test the robustness of the results 
presented in [29] , the complete expansion history of the 
background FLRW universe for i?" gravity, which resem- 
bles the ACDM model needs to be found. This is achieved 
by integrating the dynamical systems equations so that 
an orbit representing a cosmic history passes close to the 
matter dominated point G and eventually tends towards 
the late-time attractor C. In [17] it was shown that if 
1.36 < n < 1.5, G and G respectively represent a decel- 
erated matter dominated solution and late-time power- 
law acceleration. The equations for Am and TZ can then 
be solved for this background to obtain a matter power 
spectrum which can be directly compared with the one 
found for ACDM. 

This paper is organised as follows: In Section II we 
give the basic equations governing the background evo- 



lution of FLRW models and briefly present the key vari- 
ables needed to study cosmological perturbations in f{R) 
gravity using the covariant approach. In Section III we 
discuss the cosmology of i?" gravity and describe their 
main features by recasting the cosmological equations 
as an autonomous system of first order equations. We 
then integrate these equations for different values of the 
parameter n using initial conditions in the matter dom- 
inated epoch which have Hubble and deceleration pa- 
rameters equal to their ACDM values. For such initial 
conditions, we find that it is impossible to have cosmic 
histories that simultaneously have present day values of 
the Hubble and deceleration parameters close to their 
ACDM values today. In Section IV we determine the 
matter power spectra for the expansion histories given in 
V and compare them with what is obtained by integrat- 
ing the perturbation equations for the matter dominated 
solution G. We find that although the broad large and 
small scale features of the power spectrum are largely 
the same as in [30], the scale-dependent features are no 
longer present when the complete background expansion 
history is considered. We then use the observed mat- 
ter power spectrum based on the luminous red galaxies 
in the Sloan Digital Sky Survey (SDSS) [33] to directly 
constrain this class of f{R) theories of gravity and find 
that only for very special initial conditions, the models 
considered give power spectra in the Sloan wavenumber 
interval, which are in good agreement with the available 
data for luminous red galaxies. 

Finally in Section V we discuss the results and give an 
outline of future work to be done. 

Unless otherwise specified, natural units (fi, = c = 
ks = SttG = 1) will be used throughout this paper, 
Latin indices run from to 3. The symbol V represents 
the usual covariant derivative and we use the — , -f , -|- , -f 
signature and the Riemann tensor is defined by 

R\cd = W\d^c - W\cA + W\dW\e - WhcWdf ,(1) 

where the W^bd are the Christoffel symbols (i.e., sym- 
metric in the lower indices), defined by 

W\d = ^5"" {9be,d + ged,b - 9bd,e) ■ (2) 

The Ricci tensor is obtained by contracting the first and 
the third indices 
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Finally the action for /(i?)-gravity can we written in 
these units as: 
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(4) 



where R is the Ricci scalar, / is general differentiable (at 
least C^) function of the Ricci scalar and Cm corresponds 
to the matter Lagrangian. 



II. THE COSMOLOGICAL EQUATIONS FOR 
f{R) GRAVITY 

The Background FLRW equations 

In a FLRW universe, the non-trivial field equations 
lead to the following equations governing the expansion 
history of the Universe: 



3iJ + 3i/2 = __L 
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i.e., the Raychaudhuri and Friedmann equations [34]. 
Here H is the Hubble parameter, which defines the scale 
factor a{t) via the standard relation H = d/a, the Ricci 
scalar is 



R = 6H+12H^ 



(7) 



and /',/" and /'" abbreviate ^"//(Si?)" for n = 1..3 
respectively. The energy conservation equation for stan- 
dard matter 



fJ-r. 



-■iH^im (1 + w) 



(8) 



closes the system, where if is the barotropic equation of 
state. 

Note that the Raychaudhuri equation can be obtained 
from the Friedmann equation, the energy conservation 
equation and the definition of the Ricci scalar. Hence, 
any solution of the Friedmann equation automatically 
solves the Raychaudhuri equation. 



Density perturbations in f{R) gravity 

Density (scalar) perturbations may be extracted from 
any first order tensor Tab orthogonal to the 4 -velocity 
M° by using a local decomposition [28], so that repeated 
application of the operator Da = h^a^b on Tab extracts 
the scalar part of the perturbation variables. In this way 
we can define the following scalar quantities 
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-D-fi„ 



Z^Sa^D^H, C = S''D^R, 



n = a^D^R, n^a^D^R, 



(9) 



where hab is the projection tensor into the rest-spaces 
orthogonal to u°. A^, Z respectively represent the fluc- 
tuations in the matter energy density fXm and expan- 
sion Q, and TZ, SR determine the fluctuations in the Ricci 
scalar R and its momentum R. This set of variables com- 
pletely characterises the evolution of density perturba- 
tions. Then, using eigenfunctions of the spatial Laplace- 
Bcltrami operator defined in [28]: D^Q — —^Q, where 



k = 2-Ka/X is the wavenumber and Q = 0, we can expand 
every first order quantity in the above equations, so for 
example in the case of A^, we have 



A„,(i,x)=^AW(i)gW(x) 



(10) 



where ^ stands for both a summation over a discrete 
index or an integration over a continuous one. In this 



way, it is straightforward, although lengthy, to derive a 
pair of second order equations describing the fc*'' mode 
for density perturbations in f{R) gravity. They are: 
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(11) 
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Already on super-Hubble scales, k/aH <C 1, a number 
of important features can be found which allow one to 
differentiate from what is obtained in GR [29]. Firstly, 
it is clear that the evolution of density perturbations is 
determined by a fourth order differential equation rather 
than a second order one. This implies that the evolu- 
tion of the density fluctuations contains, in general, four 
modes rather that two and can give rise to a more com- 
plex evolution than the one of GR. Secondly, the pertur- 
bations are found to depend on the scale for any equa- 
tion of state for standard matter (while in GR the evolu- 
tion of the CDM perturbations are scale-invariant). This 
means that even for dust, the evolution of super-horizon 
and sub-horizon perturbations are different. Thirdly, it 
is found that the growth of large density fluctuations can 
occur also in backgrounds in which the expansion rate 
is increasing in time. This is in striking contrast with 
what one finds in GR and would lead to a time- varying 
gravitational potential, putting tight constraints on the 
ISW for these models. 

Let us now turn to the case of a general wave mode 
k. One of the most instructive ways of understanding 
the details of the evolution of density perturbations for 
a general k is to compute the matter transfer function 
r(fc), defined by the relation [35] (A,„(ki)A™(k2)) = 
T(fci)A(ki -I- k2), where k^ are two wavevectors charac- 
terizing two Fourier components of the solutions of (11) 
and T(ki) = T(fci) because of isotropy in the distribu- 
tion of the perturbations. This quantity tells us how the 
fluctuations of matter depend on the wavenumber at a 
specific time and carries information about the amplitude 
of the perturbations (but not on their spatial structure). 



In GR, the transfer function on large scales is constant, 
while on small scales it is suppressed in comparison with 
the large scales (i.e., modes which entered the horizon 
during the radiation era) [36] . In the case of pure dust in 
GR the transfer function is scale invariant. Substituting 
the details of the background, the values of the param- 
eter n, the barotropic factor w and the wavenumber k 
into (11) one is able to obtain T(fc) numerically. 

One can easily see from expressions (11) that the mat- 
ter power-spectrum in f{R) gravity theories is further 
processed after equality and would differ from the stan- 
dard ACDM power spectrum p^CDM .^^j^g^ evaluated to- 
day. The latter is widely assumed to represent accurately 
the evolution of perturbations till radiation-matter equal- 
ity, since before that the effects of any modification to the 
usual Concordance Model needs to be negligible in order 
to preserve the cosmological standard model predictions 
in the radiations-dominated epoch such as the primordial 
light elements abundances during Big Bang Nucleosyn- 
thesis. 

Therefore, these two power spectra, when evaluated 
today, would be related linearly by a transfer function 
T(fc) given by 

TjfiR) 



T^k)P^CDM 



(12) 
is obtained from the 



where T(fc) ex [A^J^^,^^ and A^ 
system of equations (11). 

On linear scales, P^ will in general depend on both 
the /(i?)-model and the scale fc, therefore differing from 
the ACDM model, where it is scale-invariant. 

In the GR limit: f{R) = R, (11) reduces to the stan- 
dard equations for the evolution density perturbations in 



GR: 
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= 0, (13) 
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(14) 



and one can easily see that the hnear evolution of CDM 
density perturbations for sub-Hubble (fc 3> aH) scales in 
ACDM is given by the well-known result: 
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(15) 



where "H = a' /a and prime (only for this equation) de- 
notes derivative with respect to conformal time. 

Notice that according to (15) the evolution of the 
Fourier modes does not depend upon k. This means 
that for ACDM on sub-Hubble scales, once the density 
contrast starts to grow after matter-radiation equality, 
evolution only changes the overall normalisation of the 
matter power-spectrum P(fc), but not its shape. 



III. DETERMINING THE EXPANSION 
HISTORY FOR i?"-GRAVITY 

To proceed, we need to fix our theory of gravity. The 
simplest and most widely studied form of f[R) gravita- 
tional theories is f{R) = aH^{R/H^)'^ , where a = a{n) 
is a non-dimensional coupling constant and Hq is the 
ACDM value of the Hubble parameter today. For this 
class of models the cosmological equations associated 
with a FLRW universe are particularly easy to analyse. 

However the aim of this investigation is to show that 
studies of different /(i?)-gravity models that share sim- 
ilar background expansion history with ACDM can in 
principle provide useful constraints on the viability of 
these models via the power spectra of matter density per- 
turbations they produce. 

The first step in the implementation of the Dynamical 
Systems (DS) approach for determining the expansion 
history for i?" gravity is the definition of the key DS 
variables. Following [17], we introduce the dimensionless 
variables: 



R{n - 1) 
HR 
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R{1 - n) 
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n,^ 
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inaH'^R" 



t(16) 



In terms of these variables, the Fricdmann equation (6) 
takes the form 



1 



y 



n. 



0. 



(17) 



An autonomous system of ordinary differential equations, 
which are equivalent to cosmological equations (5-8) can 



be obtained by differentiating (16) with respect to red- 
shift z. Here we give the equations for dust {w — 0), 
while those for a general barotropic equation of state w 
are presented in [17]: 
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(18) 



The dimensionality of the resultant system (18) can be re- 
duced further using the Fricdmann constraint (17). The 
evolution of the Hubble parameter can then be deter- 
mined by writing (3) in terms of the DS variables: 



(1 + ^) 



dh _ h{2 + ny) 
dz n — 1 



(19) 



where h = H/Hq. Furthermore the deceleration param- 
eter can be determined directly from y: 



ny 



1 



1 



(20) 



In [17] it was shown that these equations admit a number 
of fixed points of which two are particularly interesting. 
The points, labeled G and C in [17] correspond to two 
cosmologically interesting exact solutions: G corresponds 
to a matter dominated saddle point which in the case of 
dust has a ~ a^t'^"'^^ and G is the late-time attractor 



ooi 



(l-n)(2n-l) 



In [17] it was also shown that C 



with a 

and G respectively represent decelerated and accelerated 
phases of the Universe with positive energy density if n 
lies in the range 1.36 < n < 1.5. 

With this in mind let us integrate (18), by fixing the 
initial conditions for the DS variables (16) to be identi- 
cal to their ACDM values deep in the dust dominated 
era (at a redshift zq = 2000) and determine the expan- 
sion history for R" models with eight different values for 
the exponent n > 1 between n = 1.1 and 1.4 in order 
to allow the possibility of late-time acceleration. In this 
way we can determine for which values of n we obtain 
present day values for q(z) and H{z) consistent with the 
ACDM model. It is clear from the results in Table I that 



n 


1.1 


1.2 


1.27 


1.29 


1.3 


1.31 


1.33 


1.4 


ho 


0.63 


0.73 


0.94 


0.95 


1.08 


1.44 


3.27 


45.03 


qa 


0.38 


0.21 


0.15 


0.13 


0.22 


0.27 


0.21 


-0.17 



Table I: Present-day values of the Hubble ho = H (today) / Ho 
and deceleration (go) parameters for the R" models under 
consideration. Ho corresponds to the ACDM Hubble param- 
eter value today. Only n — lA provides acceleration at the 
present time, whereas n = 1.29 gives the closest value for ho 
to ACDM. 

it is not possible for i?" gravity to admit FLRW cosmic 



n exponent 


1.1 


1.2 


1.27 


1.29 


1.3 


1.31 


1.33 


1.4 


x' 


15.1394 


13.1839 


13.0184 


13.0093 


13.0104 


13.0098 


13.0102 


13.0128 


a exclusion 


1.9849 


1.4086 


1.3486 


1.3452 


1.3457 


1.3454 


1.3456 


1.3465 


% suppression 


29.5 


7.94 


4.75 


4.45 


4.37 


4.33 


4.35 


4.47 



Table II: Fits to the SDSS data for ¥C cosmology by using set of initial conditions III: eight different values of exponent n 
were investigated from n — 1.1 to 1.4. Values for x ^nd the confidence region a are presented in the second and third rows 
respectively. The data to be fitted by the theoretical spectra are taken from [33] and normalisation from WMAP7 was imposed 
for all the studied models. The fit provided by ACDM (x = 11.1996) is not improved by any of these parameter values. The 
final row gives the suppression in the overall initial amplitude required to get the best fits. For all the values, this suppression 
turns out to be smaller than 30% and is therefore in the experimental uncertainty interval for this quantity. One can see that 
the best fit corresponds to the value n = 1.29 with a suppression of 4.45% and good fits are also obtained for n = 1.3, 1.31 and 
1.33 with similar suppressions. 



histories that simultaneously have present-day values of 
the Hubble and deceleration parameters which are close 
to their ACDM values today if initial conditions are cho- 
sen in order that the expansion history is close to the 
ACDM model at early times. Even without considering 
LSS data, this already put severe constraints on these 
models. 



IV. THE MATTER POWER SPECTRUM AND 
SDSS CONSTRAINTS 



Let us now turn to the matter power spectrum. 

Taking the dominant component to be dust, the system 
(11) can be written in terms of the dynamical variables: 
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(21) 
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where prime denotes derivative with respect to redshift 
and the dimensionless quantities T?.*"' = TZ'^ /Hq and k = 
k/Ho have been introduced. Note that equations (21) 
are valid only for n 7^ 1. 

SDSS correlation data from LRG have been used to 
test the predictions from the ACDM power spectrum ob- 
tained from linear perturbation theory with WMAP cos- 
mological data to high accuracy (x^ ~ 11.2, degrees of 
freedom {d.o.f.) = 14) [33]. In what follows, we will do 
the same for this class of f{R) theories of gravity. 

To do this we first determine the cosmological back- 
ground evolution as described in the previous section 
and then use these results to solve the system of equa- 
tions (21) in order to obtain the density contrast today. 
Then, by applying expression (12) to these results, one 



can obtain the fully processed power spectra P^ for 
the above models, which can be compared to the ACDM 
predictions and the LRG data. 

Before proceeding, let us mention that three sets of 
different initial conditions were considered for the sys- 
tem (21) in order to determine how sensitive the final 
processed power spectrum is to changes in these values: 

. I: A:;Jo = Te^lo = 10-5, A:;;|o = 7^^•'|o = IQ-^, 
. II: Af;jo ^Te'^lo = lo-^ Aj;;|o =7^'^'|o = 10-8, 
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where the subscript refers to the initial redshift zq 
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Figure 1: The left panels show the transfer function T{k) = |Afe/A^'-^^^(2: = 2000) p evaluated today {z = 0) for wavenumber k (in 
hMpc~^ units) in the range 0.0001 to 0.1 for the initial conditions sets I, II and III as described in the bulk of this investigation. The 
transfer functions T(fe) on left panels have been normalised in such a way that the curves coincide on large scales. One can see how for 
conditions I and II the power spectra are hardly distinguishable whereas for the conditions III the results differ. On the right panels 
we present the corresponding linear matter power-spectra P{k) for ACDM and R" models for n = 1.1, 1.2, 1.27, 1.29, 1.3, 1.31, 1.33 
and 1.4 data correspond to SDSS data [33]. All the power spectra were assumed to have an arbitrary overall normalisation at the scale 
k = 0.01/iMpc~^ in order to find the best fit to the data. Conditions I and II lead to power spectra in complete disagreement with the 
observed data. Conditions III, due to the almost fiatness of the spectra in the range covered by Sloan data present a good fit to the data. 
On the right bottom panel we show in a window the relative discrepancy between the ACDM and the i?" fits power-spectra for every 
studied exponent. The smallest discrepancy in scales A: > 3 ■ 10~^ happens for n = 1.1 whereas for smaller scales, all the remaining values 
of n provide similar relative error around 5 ■ 10"'^. 



the system (11) can be understood as providing scale- 
invariant initial conditions for the variables A^ and TZ^ 
and their first derivatives are all taken to be small (but 
non-zero) at the initial redshift. On the other hand, set 
III we set the first derivatives of A^ and IV' to zero at 
z = zq. This choice has important consequences for the 
obtained spectra. 

For each value of n, we present in Fig. 1 both the 
transfer functions and the processed power spectra for 
all the initial conditions sets just mentioned. 

One can see how for these models that the transfer 
functions have a nearly-flat plateau on large scales [30] 
regardless of the set of initial conditions. On intermediate 
scales however, the density contrast behaviour and its 
amplitudes today depend both upon the value of n and 
the initial conditions. The left panels in Fig. 1 clearly 
illustrate this. 

Then, by using expression (12) for the transfer func- 
tion, one can then obtain the processed power spectra 
and compare these theoretical results with the LRG data. 
We find that initial conditions I and II are not able to 
provide good fits to the SDSS data due to the fact that 
in the required scales (0.01 — 0.1 h Mpc~ ) the spectra 
are not flat but decrease with the wavenumber k (see the 
left panels in Fig. 1). On the other hand, initial condi- 
tions III give power spectra which are in good agreement 
with the LRG data. This is due to the almost-flat trans- 
fer function in the LRG data range. Note however that 
the initial amplitude was assumed to be a free parameter 
which was determined to achieve the best flt. 

For this set of initial conditions. Table II give the x^ 
{d.o.f. = 14) analysis for the eight studied i?" models 
when their respective spectra evolutions are fitted to the 
SDSS data, assumed to be non-correlated. We also in- 
clude the value for the confidence regions a with respect 
to ACDM as well as the overall amplitude suppression 
in the initial scale (k = 0.01 h Mpc~ ) to get the best 
fits after a Least Square method analysis. It is clear that 
none of the i?" models under consideration acquire the 
same goodness of fit as the ACDM model as seen in Table 
II in the a exclusion regions. 

For completeness and in order to emphasise the im- 
portance of using the complete expansion history for the 
background, we have also given |A^p for models whose 
background evolution is given by the exact saddle point 
solution G in the case of dust {w = 0). We do this 
for the same parameter values n and initial conditions 
I, II and III as shown in Fig. 2. These results agree 
with previous investigations [30] that showed how when 
this background scale factor is assumed the spectrum is 
composed of three parts corresponding to three different 
evolution regimes for the perturbations. In this scenario, 
on intermediate scales the interaction between the two 
fluids (dust and curvature) is maximised and the curva- 
ture fluid acts as a relativistic component whose pressure 
is responsible for the oscillations and the dissipation of 



the small scale perturbations in the same way in which 
the photons operate in a baryon-photon system [30]. If 
we compare these results to the left panel in Fig. 1, 
we conclude that the scale-dependent features in Fig. 2 
are washed out when the complete background expansion 
history is considered, however the main large and small 
scale features of the power spectrum found in [30] are 
retained. 



V. DISCUSSION AND FUTURE WORK 

In this paper we presented a complete analysis of the 
background and matter perturbations for one of the most 
widely studied modified gravity theories: i?" gravity with 
n > 1. Both the cosmological background evolution and 
linear perturbation equations were solved by combining 
the dynamical systems approach for the background and 
using the 1-1-3 covariant approach to evolve the mat- 
ter perturbations, without assuming any intermediate 
(quasi-static) approximation. 

We solved the background equations for different val- 
ues of the parameter n using initial conditions in the 
matter dominated epoch, with Hubble and deceleration 
parameters equal to their ACDM values. For such initial 
conditions, we found that it is impossible to have cos- 
mic histories that simultaneously have present day values 
of theses cosmological parameters close to their ACDM 
values today. In fact, of the eight models considered, 
only n ~ \A provided a negative deceleration parameter 
today, but gave a Hubble parameter completely incom- 
patible with its observed value, while n = 1.29 gave the 
closest value for the present-day Hubble parameter to 
ACDM but exhibits no late time acceleration. 

We then used the observed matter power spectrum 
based on the luminous red galaxies in the Sloan Digi- 
tal Sky Survey to further constrain these models. For 
the studied exponents, we found that all the models 
gave rise to an almost-flat transfer functions in the Sloan 
wavenumber interval provided very special initial condi- 
tions are chosen. In this case the best flt to the data was 
found for the value n = 1.29 with a suppression of 4.45% 
and good flts were also obtained for n = 1.3, 1.31 and 
1.33. The exponent n = 1.4 (the only one providing ac- 
celeration today) required a suppression slightly bigger 
(4.47%). None of the studied exponents were however 
able to fit the data as well as the ACDM model. 

It is clear from this analysis that i?" gravity does not 
successfully meet any of the cosmology requirements for 
it to be considered as a viable alternative to the standard 
model. This work does however illustrate the utility of 
our approach and it should be possible to use these tech- 
niques to constrain which f(R) theories remain consistent 
with current data even if they are indistinguishable from 
the ACDM model at the level of the FLRW background. 
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